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PROCESS FOR DETERMINING A THR^S^ 
A TWO DIMENSIONAL IMAGE, PARTICULARLY THE STRUCTURE OF A 

BONE 

5 The present invention relates to a process for 

determining a three dimensional structure from a two 
dimensional image, particularly that of a bone. 

Osteoporosis is a malady which affects a large number 
particularly of women after menopause toward age 50, given 

10 that this malady can affect anyone at any age. It is 
characterized by a low bone mass and a deterioration of the 
bone tissue. These degradations lead to great risk of 
fractures, particularly of the hip, of the vertebral column 
and of the wrist. 

15 of course there exist a number of risk factors but 

this is not enough to determine the probability of a person 
fracturing an element of his skeleton given that there is 
no preliminary symptom. 

This malady is widespread in the population and will 

20 increase because of particularly the sedentary lifestyle 
and the aging of the population. 

However, the treatment of persons after fracture is 
very long and very costly because this involves long term 
care. The consequences are often serious because they lead 

25 to invalidism or even death of the patient. 

Also, it is necessary to prevent this type of 
affliction but after determining the fracture risk, which 
is the only resort of the practitioner. 

If it is possible to predict the risks, the 

30 practitioner then has treatments for hormonal therapy and 
by other pharmacological products such as calcitonin or 



biophosphenates, in addition to advice as to healthy 
lifestyle to attempt to avoid this affliction. 

There are apparatuses to determine the bone density, 
called densitometers. 
5 It has been thought that if the bone density were 

known, and compared with standard curves, the practitioner 
could vthen determine the risk and establish a suitable 
diagnosis and prescribe a satisfactory treatment. 

This is what is done at present. 
10 From U.S. Patent 5, 774, 520, it is known that the 

probability of a fracture is directly connected to bone 
density. 

Densitometers available commercially use x-rays or low 
energy gamma rays. The absorption being proportional to 

15 the bone density, the two can thus be correlated to have a 
satisfactory result permitting comparisons. This patent 
teaches working with two sources of the emission of photons 
so as to process the dispersion that arises because of the 
osseous materials being considered as homogeneous, but the 

20 other liquid materials, the muscular tissue, the cartilages 
cannot be considered as a single homogeneous material. 

Another U.S. Patent 6,385,283 uses the density but 
combines this measurement with images permitting 
determining the risk of fracture. These images are made 

25 from the vertebral column of the patient and the operator 
determines the presence of beginning fractures. Moreover, 
as in the known prior art, this information is compared 
with data acquired from numerous cases. This diagnosis is 
refined by the practitioner from the history of the patient 

30 and by introduction of supplemental risk factors. 

Patent application WO 86/07531 proposes providing an 
image of a given bone, for example the calcaneum, at an age 
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at which the patient normally has his full osseous mass. 
Then, the process consists in providing successive images 
of this same bone to permit the practitioner to carry out 
comparisons and to determine the development so as to draw 
5 conclusions as to the degree of risk of fracture. 

All these methods and other devices are based on a 
single measurement of the density of the osseous material. 

However, so that the practitioner can arrive at a 
satisfactory diagnosis, it is necessary that he have 
10 information other than the density alone. 

The definition of osteoporosis is given by the WHO: 
"disorder characterized by a low osseous mass and 
alterations of the microarchitecture of the osseous tissue, 
leading to increase of fragility of the bone and as a 
15 result an increase in risk of fracture" 

Thus there exist analyses of osseous material in vivo 
that are more sophisticated, from analyses of blood and 
urine useful as osseous markers but these analyses are 
costly, difficult for the patient and, in any case, give 
20 indications of composition of the material but not 
indications of its structure. 

However, it will be understood that the osseous 
material has a different resistance according to whether 
there is more or less quantity of material, which results 
25 from the density measurement. On the other hand, for a 
same osseous density and hence for a same quantity of 
material, according to the architecture of the osseous 
structure, the mechanical resistance can vary greatly. 

What are useful for a practitioner, are the mechanical 
30 properties of the bone of a patient leading directly to the 
risk of fracture. When there is a decrease of osseous 
resistance, there is a conjugation of two factors, the 
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decrease of the osseous mineral quantity, and the change of 
the bone structure. 

At present, there is no means permitting determining 
the bone structure. There can be used magnetic resonance 
5 imagery or scanning, but these are difficult and costly 
examinations, particularly if it is necessary to carry out 
a longitudinal study of the patient. 

Moreover, the capture of three dimensional images from 
imaging means is at present limited by the capacities of 
10 these means because the spatial resolution is greater than 
the dimensions of most of the connection spans permitting 
carrying out connectivity, which is to say the connections 
between the three dimensional network nodes. 

The process according to the present invention will 
15 now be described in detail, so as to permit by analysis of 
a two dimensional image, estimating the mechanical 
parameters . 

According to the invention, the process of 
determination of the mechanical resistance of a bone, 

20 according to the invention, from a digitized two 
dimensional image, obtained by imaging, is characterized in 
that there is carried out a correlation between the bone 
mineral density determined from this two dimensional image 
by any means suitable for this type of imaging, and a 

25 structural parameter obtained from the same two dimensional 
image . 

Particularly, recourse is had to a correlation of the 
exponential type. 

This correlation is used by associating the bone 
30 mineral density and said structural parameter, to determine 
the ultimate strain C u of the bone. 
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More precisely, the structural parameter a is 
determined, obtained by the following sequence of steps: 

a) choosing at random a pixel of the two dimensional 
image which is at a gray level h(0), 
5 b) choosing a straight line from this point having a 

direction also determined at random, 

c) moving a distance r along this straight line, h(r) 
being the gray level of this new point, 

d) computing the variance of the gray levels with the 
10 formula: V(r) = [h(r) -h(O)] 2 , 

e) tracing the curve associated with V(r) on a log-log 
scale, and 

f) determining the slope of this log-log curve which 
represents said parameter a. 

15 So as to improve the precision, steps a) to d) are 

repeated a sufficiently large number of times to make the 
function variance V(r) converge, by means of the assembly 
of repetitions. 

According to another characteristic of the invention, 
20 there is carried out a correlation between the bone mineral 
density obtained from this two dimensional image and said 
parameter a evaluated from the same two dimensional image 
according to the mathematical model: 

C u ' = b 0 + bi* exp (b 2 * DMO) * a 
25 - wherein b 0 , bi, b 2 are the coefficients obtained 

by nonlinear regression and C u ' the prediction 
of the ultimate stress C u of the bone. 
There is determined a correlation between the 
parameter a and a three dimensional parameter of the 
30 trabecular network and the bone and a three dimensional 
parameter can be the connectivity density Xv- 
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To support this description, drawings are attached and 
the figures that they represent illustrate the description 
and show essentially the results obtained. 

Figure 1: curve of results with a linear model, C u as 
5 a function of C u ' 

Figure 2: curve of results with an exponential model, 
C u as a function of C u ' 

Figure 3: graphical representation of C u as a function 

of Xv 

10 Figure 4 : representative of f s as a function of Xv 

Figure 5: curve of the function V(r) 

Figure 6: determination of the parameter a from the 
log-log curve of the function V(r) 

Figure 7: bone mineral density curve as a function of 

15 f s 

Figure 8: curve of a as a function of Xv 

Figure 9: curve of bone mineral density as a function 

of a 

Figure 10: linear model C u as a function of bone 
20 mineral density 

Figure 11: curve of C u as a function of C u ' . 
The study related to a trabecular or spongeous bone 
which constitutes about 20% of the bone material. The 
cortical bone which surrounds it ensures the rigidity of 
25 the entire bone and hence of the skeleton, whilst the 
tribecular bone ensures the absorption and resistance to 
compressive forces . 

According to the invention, it is considered that the 
tribecular structure comprises spans which are 
30 interconnected and the mechanical resistance of the whole 
results from this connection, which is to say the number of 
closed loops. The mechanical resistance also results from 
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the degree of mineralization of these spans, mineralization 
which is an important parameter. 

The process consists in analyzing a digitized image 
obtained in the present case from X ray emission. This 
5 image is a projection on the surface of a three dimensional 
structure having been traversed by the emitted photons. 

Each pixel which constitutes the image can be analyzed 
independently but then there can be determined only the 
bone mineral density although it would be desirable not 
10 only to analyze a pixel in an isolated fashion but to 
analyze them each with respect to the others. 

The radiation spectrum should be the most 
monochromatic possible so as to avoid too great drift of 
the energy of the photons. 
15 Thus, if the energy varies the penetration will vary 

and will give nuances of gray on the image that can lead to 
errors by giving the impression of absorption there where 
there is none. 

From this image obtained with a narrow spectrum, it is 
20 known that each elemental volume undergoes a direct action 
of the photons, perpendicularly but also an indirect action 
of the incident rays diffused by the adjacent elemental 
volumes after they have been themselves traversed. 

So as to overcome diffusion, there can be used 
25 preferably two separate energies. There are two equations 
with two unknowns and there can thus be deduced the portion 
of real attenuation which directly interests the present 
process and the portion resulting from diffusion. 

If the ultimate stress C u (MPa) is sought, it is 
30 necessary to take account of the following parameters: 

V 0 : analysis volume (mm 3 ) 
V s : solid volume (mm 3 ) 
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f s : volume fraction of solid 
with the correlation f s = V s / V 0 

(3 0 : number of mass of the solid portion (which 
is to say the assembly of the connected or 
5 disconnected portions) , this number is 

generally equal to 1 . 

p 2 : number of internal surfaces, which is to 
say the holes resulting from the internal 
microporosity of the spans or the working 

10 scale, this porosity not being visible, this 

number is generally equal to 0. 
Nep : Euler-Poincarre number 
X : connectivity 
with the relationship x = Po + P2 - N E p 

15 - Xv: connectivity density (mm" 3 ) 

with the relationship Xv = X / v o 

The process according to the present invention seeks 
to permit correlating two measurable parameters and to find 
the relationship which connects them. 
20 1/ If there is used a linear model such as: 

C u ' = a 0 * f s + ai 
C u ' being a predicted value for C u 
a 0 and ai are coefficients of linear regression, 
there is a simple correlation of C u with f s . 
25 The obtained curve is shown in Figure 1 and it will be 

seen that the dispersion is very great. If for example f s 
is taken to be constant, the variation of the values of C u 
is great. 

This linear model cannot be followed, because its 
30 results are insufficiently precise, but it can nevertheless 
serve as a comparison. 

2/ If there is used an exponential model such as: 
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C u ' = b 0 + bi * exp(b 2 * f s ) * Xv 
Wherein b 0/ b x , b 2 are coefficients obtained by 
nonlinear regression (method of generalized least squares) 
and C u ' the prediction of the ultimate stress C u of the 
5 bone. 

The dispersion is greatly limited, see Figure 2. C u 
is estimated from f s and Xv 

These variations are shown in Figure 3. 

If f s and Xv increase, C u increases. This is perfectly 
10 natural, if the volumetric fraction of solid increases and 
the connectivity increases, and hence the ultimate stress 
increases . 

If Xv is fixed, then C u is a function directly of f s/ 
in an exponential manner. At constant connectivity, the 
15 ultimate stress increases as a function of the increase of 
the bone density. 

If f s is fixed, then C u is a direct function of Xv, in 
a linear manner but the ultimate stress C u decreases when 
the connectivity increases. 
20 If the solid fraction is constant, it will be seen 

that the ultimate resistance C u decreases linearly as a 
function of the connectivity. 

The more nodes in the structure, the more this 
structure becomes fragile, in a surprising manner, contrary 
25 to a well established prejudice. 

There can thus be traced the model of representation 
of Figure 4, which shows the values of C u . 

To determine connectivity, it is thus necessary to 
find a parameter which follows a same law and which will be 
30 as independent as possible of the bone density so as not to 
be influenced. 
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There could be carried out a computation from two 
measurable parameters such as the bone volume and the 
connectivity, these two parameters being obtained 
particularly by MRI . 
5 On the other hand, obtaining values by this method 

uses costly apparatus, hardly available, which leads to 
predictions of a high price preventing a longitudinal 
follow up and hence a massive expenditure by the patients. 

It is thus necessary to find two values measurable for 

10 example from simple densitometry, which is itself 
altogether approachable and which can be regularly 
repeated. Thus, densitometries have already been carried 
out in a massive way for numerous applications. The cost 
for patients can be very great. 

15 The process according to the present invention 

consists in determining these two measurable parameters as 
the bone volume fraction and connectivity, or more exactly 
parameters which could be correlated with them in a 
sufficiently narrow way to give satisfactory results. 

20 The first rapidly measurable and reliable parameter is 

the bone mineral density. Apparatus has been developed to 
carry out these measurements in a reliable and 
reproduceable way. If the bone mineral density curve is 
traced as a function of f s , it will be seen that the 

25 dispersion is low and that bone mineral density can be used 
in lieu of and in place of f s . See Figure 7. The first 
parameter is obtained. 

The second parameter is called a and is determined 
from a digitized two dimensional image for example obtained 

30 with x ray. 
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To determine this new parameter a, there is first of 
all chosen at random a two dimensional pixel of the image 
which is at a gray level h(0). 

There is then selected a straight line from this point 
5 having a direction also determined at random. 

We move by a distance r along this straight line, h(r) 
being the gray level of this new point. 

The variance of the gray levels is then calculated, 
which gives : 

10 V(r) = [h(r) - h(0) ] 2 

The process is repetitive, with a repetition of these 
steps a number of times sufficiently great, and at random, 
to determine a. There is thus caused to converge the mean 
variance function V(r) from the assembly of repetitions. 

15 This permits obtaining by computation an assembly of points 
and the associated curve is shown in Figure 5. 

Thus, the more one determines points having a same 
gray level adjacent the selected pixel, the more one finds 
himself in the material and hence the more one is 

20 homogeneous and the more can be found the connection modes. 

If the log-log curve of this function is traced, 
Figure 6, it will be seen on the first five points in this 
instance, a linear behavior and the slope of this straight 
line is the coefficient called a and suitable for the 

25 present invention. 

If there is then traced the curve a as a function of 
the volumetric connectivity Xv# Figure 8, it will be seen 
that there is a sufficient correlation so that Xv can be 
replaced by a. 

30 The curve of Figure 9 can be deduced from the 

preceding curves, as a diagram which shows the relation 
between bone mineral density, a and C u . 
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Bone marrow density and a are two parameters directly 
measurable and/or at least calculable from a same two 
dimensional image . 

It will be seen on this diagram that when the bone 
5 marrow density remains constant, the ultimate stress 
diminishes when a increases. 

This is entirely equivalent to the exponential model 
indicated above, with the bone fraction as a function of 
conductivity, see Figure 4, except that the connectivity Xv 
10 was not directly measurable whilst a is calculated from the 
two dimensional image, and hence quantifiable. 

We can then write the same relationship but with the 
parameters determinable from the same image: 

C u ' = b 0 + bi* exp(b 2 * DMO) * a 
15 This is what Figure 9 shows. The curve of Figure 11 

characterizes this exponential model of C u ' as a function of 
a. There will be seen a slight dispersion compared to that 
of the linear model whose curve is shown in Figure 10. 
This is a non-limiting example which shows the possibility 
20 of connecting the bone density whose value is determined 
from a two dimensional image, to a structural parameter 
determined from this same two dimensional image to obtain 
the ultimate stress C u . 

It will be seen that it is possible to determine a 
25 correlation between the parameter a evaluated from the two 
dimensional image and a three dimensional characteristic of 
the trabecular network of the bone as the density of 
connectivity Xv 

There can thus be given an example obtained from 
30 femoral bone images obtained by absorption measurement with 
double energy which replaces an image obtained by x rays or 
gamma rays of low energy. 
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These images are those of a young woman of 17 years 
and the other of a woman aged 69. The two dimensional 
images thus obtained are analyzed according to the process 
of the invention. 
5 There will be seen from measurement of bone mineral 

density the same values obtained directly by the 
densinometric system, in the two cases: 0.730 g/cm 2 in the 
trochanterial zone . 

The calculated parameters a are respectively 0.56 for 
10 the young subject and 0.71 for the aged subject, thereby 
showing a notable difference and the parameter a which is 
the least corresponds to the best mechanical resistance of 
the bone C u , that of the younger subject. 
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